Novel antiferromagnetic quantum phase transition in underdoped cuprates 
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We investigate a zero-temperature itinerant antiferromagnetic transition where the fermions possess 
a d-wave gap. This problem pertains to both the nodal liquid insulating phase and the d-wave 
superconducting phase of the underdoped cuprates. We find that a non-trivial quantum phase 
transition exists, and that the quantum critical point is dominated by a long-ranged interaction 
(\x — y\~ 2d ) of the Neel order parameter, which is induced by the Dirac-like fermions near gap 
nodes. We formulate a Ginzburg-Landau functional and estimate the critical exponents via the 
large-n expansion method. 
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Recent experiments have revealed that underdoped 
cuprate materials exhibit an exotic pseudogap phase un- 
der a characteristic temperature T* . There is ample evi- 
dence of the pseudogap phase above the superconducting 
critical temperature [0 ; even a low-temperature pseudo- 
gap may exist in vortex cores in the mixed states 0. 
Although the underlying microscopic mechanism of the 
pseudogap phenomena is not understood, we may never- 
theless study the effective theory and its phenomenolo- 
gies. Based on the smallness of the superfluid density 
in underdoped cuprates, the pseudogap phase may be 
described as a superconductor whose phase coherence 
is destroyed by strong thermal phase fluctuations but 
where the gap amplitude is robust Hf|. As a possi- 
ble zero-temperature analog of the pseudogap state, the 
idea of nodal liquid has been proposed, as a supercon- 
ductor which is quantum disordered by quantum phase 
fluctuations [||. The nodal liquid is a non-magnetic in- 
sulating spin-liquid state, which can be effectively de- 
scribed by charge-neutral fermionic quasiparticles with a 
d-wave gap, but with a non-zero insulating gap in terms 
of the charge-ful electrons. This exotic state can be an in- 
tervening phase between the undoped antiferromagnetic 
(AF) insulating state and the the d-wave superconduct- 
ing state. Further experimental and theoretical studies 
on the insulating and AF phases are needed to confirm 
this idea. 

In this report, a distinctive consequence of the nodal 
liquid state is presented: We show that quasiparticles 
with a d-wave gap cause a characteristic critical behavior 
of quantum phase transitions in the nodal liquid phase, 
which is distinct from that of the conventional Fermi liq- 
uid More specifically, we study an itinerant AF 
transition at zero temperature. We expect, however, that 
without much modification the same is true of certain 
classes of itinerant magnetic transitions in a supercon- 
ducting state of unconventional symmetry. The origin of 
the characteristic scaling is the long-ranged interaction 
between the order parameter which is induced by mass- 
less Dirac-like fermions near gap nodes; after integrating 
out the fermionic degrees of freedom, a term of the form 
N(x) • N(y)/\x — y\ 2d is induced, where N is the Neel 



order parameter, d is the spatial dimension, and x, y are 
(d + l)-dimensional space-time vectors. We show that 
a non-trivial quantum critical point exists at the transi- 
tion and that the abovementioned long-ranged interac- 
tion dominates the critical behavior J§|. The same argu- 
ment holds without modification for a zero-temperature 
AF transition in the d-wave superconductor although it 
is uncertain whether such a transition takes place in the 
superconducting state of cuprates. This long-rangedness 
of the order parameter is similar to the case of itinerant 
ferromagnetic transition j^] where non-critical soft modes 
of the Fermi liquid induce such non-trivial interactions. 
The difference is that in the case of Dirac fermions, ad- 
ditional soft modes are not needed to generate the long- 
ranged order parameter interactions. 




FIG. 1. Schematic geometry of the Fermi surface and the 
d-wave gap. The solid line is the nested Fermi surface, and the 
dashed line represents the gap in the quasiparticle spectrum. 
Q depicts an AF ordering wave- vector. 

First, we construct an effective theory of nodal 
fermions with commensurate AF exchange interaction. 
We will not discuss AF ordered phase in this paper. 
Here and throughout we do not make a distinction be- 
tween quasiparticles in the d-wave superconducting state 
and charge-neutral nodal quasiparticles in the nodal liq- 
uid phase. We consider a system close to half-filling 
where the gap nodes are located on the nested Fermi 
surface. (See Fig. Q). Since we are only interested in 
the low-energy properties, we integrate out the fast mo- 
mentum degrees of freedom and retain only the momenta 



1 



near the Fermi surface. Assuming that we can construct 
such a low-energy effective theory of the Hubbard model 
near half-filling |10|, we may express the effective action 
S = Sq + Saf in the following form: 



ui k 

+A^ T ( k ;^) ip L (-k; -gj) + c.c] , 
Saf = ^J |fr^n Q (q,T)-iiQ(-q,T), 



(1) 



(2) 



where n Q (q,r) = E k 3 ^( k > ^^^(k + Qk + q, r), 
and Ak is a d-wave gap. Here |q| <C kp and Qk is a 
commensurate nesting vector antiparallel to the Fermi 
velocity at k. Also we take J < so that an AF transition 
takes place. 

Before we discuss the AF critical point, we revisit re- 
sults of the mean-field theory. To obtain a mean-field 
theory, it is sufficient to employ a mixed representation 
of the momentum to separate the Fermi surface label and 
the small relative momentum near the Fermi surface as 
follows: 



S = So 



d x dr 



- — N 2 (x) 
2J V ' 



(3) 



- £ e-*V x N(z) • i^t(p; x) 5^(p + Q p ; x) 



where 



u p,q 

+a;v t (p;?) ^ i (-p;-9) + cc.} . 

Here we have separated the AF exchange interaction by 
introducing a staggered magnetization N via Hubbard- 
Stratonovich transformation. The momentum p labels 
points on the Fermi surface, and q is a small momentum 
deviation from the Fermi surface, q — (q, uS) is a (d + 1)- 
dimensional wave- vector and similarly, x = (x, t) . In this 
representation, the d-wave symmetry of the gap dictates 
that A r 



A_p and A p 



-A p+Qp . Taking |N| 



iV z =constant, we obtain, for instance, the quasiparticle 
Greens function Q{p;q) = -(iu; + £ q )/(w 2 + £ 2 + | A p | 2 + 
N 2 /A). As expected, the quasiparticle gains a gap of the 
form ^/|A p | 2 + iV 2 /4. 

From the following mean-field self-consistency condi- 
tion, 

W\=T,\ {^TMi^x) vAp + Q p ;*)> , (5) 



we obtain the relation 



\N\/2- 



A 2 , 



7V 2 /4 = Eac- 



-i/\J\n 



(6) 



where Ea is an upper energy cutoff which is of order 
Ef, N is the density of states at the Fermi surface, and 
A is the maximum gap. This suggests that due to the 
gap in the quasiparticle spectrum, the AF transition is 
inhibited to such a degree that a sufficiently strong ex- 
change interaction is needed for a transition. In fact, 
we infer from Eq. (||) that the AF coupling J needs to 
be at least of comparable strength to the superconduct- 
ing pairing strength; on a circular Fermi surface, even 
stronger exchange interactions would be required ]Tl| . 
This is in contrast to the case of Fermi liquid on a nested 
Fermi surface where arbitrary strength of exchange inter- 
action leads to an AF instability fllO|| . Therefore, a weak- 
coupling fermionic renormalization group (RG) analysis 
is inadequate for this problem. Instead, we turn to the 
Ginzburg-Landau (GL) functional of the Neel order pa- 
rameter below. 

Now we come back to Eq. (|l|) which reduces to a spa- 
tially anisotropic massless Dirac fermionic action near the 
gap nodes. In terms of the small momentum deviation 
q from the Fermi wave- vector, the fermion Greens func- 
tions have poles at u> « ±i^J vpq\\x 2 + vj^q^, where vf is 
the Fermi velocity, va is the slope of the gap at the node 
in the momentum space, and (q±) is the component of 
momentum parallel (perpendicular) to the Fermi velocity 
at the node. Although the bare theory does not have an 
exact Lorentz symmetry due to this anisotropy, we as- 
sume that even if we enforce Lorentz symmetry by hand, 
the quality of the quantum critical point is not crucially 
changed || and take Vf = 1>A = 1 for simplicity. 

To arrive at the GL functional of the order parame- 
ter, we integrate out the fermionic degrees of freedom. 
For now we assume that this procedure is under control, 
although we will show below that at T = it is beset 
with singularities in GL coefficients. We confine our dis- 
cussion to an AF transition sufficiently far away from a 
superconductor-insulator transition, and retain only the 
Neel order parameter. After formally integrating out the 
fermionic fields, 

5[N] = / d d xi dri...d d x n dr n T n {x\, x n ) 

xN( Xl )...N(x n ) , (7) 

where n is an even integer. The Fourier transform of T2 
JPH is obtained from 

TM ' V) = ~h + \ J S^> (k ' W) S(k + Qk + q,^)] , 

k 

(8) 

where Q is the diagonal component of the quasiparticle 
Greens function. Collecting the lowest powers of the mo- 
mentum and frequency, we find that for 2 < d < 3, which 
is the range of dimensionality of our interest, the leading 
term is non-analytic in \q\: 
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r 2 = t + c\q\ G 



(9) 



where \q\ — \/ Q 2 + v 2 and c > 0. Similarly, one can show 
that r„ ~ |(7| d_Tl+1 . Using the Gaussian approximation, 
(N-N) = l/(t + c 2 \q\ d ~ 1 )- Taking t as the distance from 
the critical point, we obtain various critical exponents: 
7 = 1, z = 1 , v = l/(d — 1), and n = 3 — d. By 
including higher order GL expansions and introducing 
an 0(3) symmetry-breaking term h • N, we can obtain 
the equation of state. In the presence of the staggered 
magnetic field, a non-zero N value will cut off the infrared 
singularities in T n by providing a low-momentum cutoff 
H . Therefore, T n ~ ]\jd-n+i obtains, and the equation 
of state (8S/5N = 0) is tN + N d w h. From this, we 
obtain two more exponents: = l/{d — 1) and 5 = d. 

Fixing the Gaussian term as marginal, we may esti- 
mate the scaling dimensions of higher GL expansions. By 
transforming q — > bq and N — > b~ d N, we find that the n- 
th GL expansion in TV scales as b n ~ 1 ~ d . This means that 
if we have a local GL theory, all high order expansions 
in N are irrelevant in the tree level. However, since the 
coefficient T n scales as b d ~ n ~ x due to the long-ranged 
interaction, all higher order terms gain marginal scal- 
ing. Therefore, there arc infinitely many marginal opera- 
tors in this theory, and the Gaussian approximation may 
break down. This is because of the singularities in r„ 
that occurs as T — ► 0, which suggests that a more careful 
RG analysis is called for. The result of the e = 4— (d+ 1) 
expansion in this model shows that there exists a non- 
trivial stable fixed point at the AF transition || . In this 
paper we use the large-n expansion method by introduc- 
ing n fermion species to estimate the leading correction 
to the critical exponents and take n = 1 at the end. 

We come back to Eqs. (0) and (||), and introduce fic- 
titious n fermion flavors which couple to the Neel order 
parameter symmetrically. We take the upper momentum 
cutoff A < Ao and re-express the effective action in terms 
of the fermions coupled to the Neel field as follows 

n ~ 

S cS =J2 S 0+ I duJ A[N(<7) ■ N(-q) (t + q 2 )} (10) 

i=i 

-g ^ e -^- x N(x) • -V^-Or) ff$^{x)} , 

j 

where each fermion species is labeled by the index j. 
We have not included higher order local expansions in N 
since more significant non-local expansions will be per- 
turbatively generated by the long- wavelength fermions. 
Here each Sq is of the same form as Eq. (|l|) except that 
the momentum and frequency are confined to within a 
sphere of radius A <C Ao so that the fcrmionic spectrum 
can be well-approximated by the Dirac spectrum. 

Now we take the large-n limit and estimate the criti- 
cal exponents by calculating the Neel propagator. Note 
that n is not the number of order parameter compo- 
nents but the number of fermion species. We find that 



this is a convenient choice since the singular GL coeffi- 
cients in Eq. (Q) are ultimately determined by the long- 
wavelength fermions and the coupling g in Eq. (|l0|). In 
this way, we can systematically calculate ^-independent 
critical exponents as expansions in l/n by the pertur- 
bation technique |l2|. In Eq. (^), each vertex function 
r„ consists of one fermionic loop which has a factor of 
71. This means that l/n expansion is equivalent to loop 
expansions in the Neel field — l/n takes the role of Ti. 
In practice, it is more convenient to directly work with 
fermions rather than to integrate them out at the outset 
to perform loop expansions in the order parameter field. 
Figure || (a) shows the leading Neel propagator in the 
large-n limit, which is (N(q)N(—q)) = [ng 2 \q\/lQ + t}^ 1 
in d = 2, and ~ + t]^ 1 in general dimensions. 

Note that the leading non-analyticity is already obtained 
here. The critical exponents at the leading order agree 
with the Gaussian approximation of Eq. (|t]). In order 
to estimate the l/n correction to the critical exponents, 
we calculate the diagrams shown in Fig. || (b). Here we 
focus on two spatial dimensions, as it is difficult to carry 
out analytic calculations in general dimensions. Setting 
t = 0, we can calculate the correction to r\ by finding the 
coefficient of the term \q \ In |g| in the Neel self-energy cor- 
rection. In d = 2, since there is fermion self-energy cor- 
rection ~ (l/n) fcln |fc|, we may expect a Neel self-energy 
correction of the form (l/n) |q|ln|g| on the dimensional 
ground. Explicit calculation shows that there is no such 
logarithmic correction from any of the three diagrams in 
Fig. I (b). Therefore, we obtain n = 1 + 0(l/n 2 ). Set- 
ting a non-zero t > 0, we may obtain the correction to 7 
by finding the coefficient of the self-energy correction of 
the form tint, and we find that 7 « l+0.78/n+O(l/n 2 ). 



WV = + 
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FIG. 2. (a) The thick wavy line denotes the leading order 
Neel propagator in the large-n limit, and the thin solid line is 
the bare fermion propagator. The dashed line is the bare AF 
exchange coupling and the bubble is the bare fermionic AF 
spin susceptibility, (b) Corrections to the Neel self-energy of 
order 0(l/n). 

The long-ranged interaction discussed above is domi- 
nant only in the T — > limit, but the finite-temperature 
crossover behavior 0,|l^] is beyond the scope of discus- 
sion in this report. We merely comment on the classi- 
cal finite-temperature transition for 2 < d < 3 where 
the finite-temperature Neel order is assumed to be stabi- 
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lized by weak interlayer coupling of cuprates, assuming 
that the nodal quasiparticles survive at non-zero temper- 
atures. When the temperature is the largest energy scale 
of the problem, the non-local interaction of the Neel order 
parameter gives subleading power laws near the transi- 
tion. As a result, the leading Gaussian GL expansion is of 
the form / d d g(q 2 +i)N(q)-N(-q) where t = (T-T c )/T c . 
This coincides with the usual classical GL functional, and 
therefore we find no distinctive critical behavior in this 
case. 

As pointed out by Hertz M, AF transitions on nested 
Fermi surfaces and superconducting transitions in Fermi 
liquids are two other examples which suffer from the 
T — > singularities in GL coefficients. In fact, in these 
two examples, zero-temperature GL expansions are more 
difficult because of In |q| singularity in the Gaussian term. 
This is related to the fact that the Fermi liquid ground 
state is always unstable in these cases. 

We have shown that the existence of a d-wave gap in 
the quasiparticle spectrum leads to distinctive critical be- 
havior at zero-temperature phase transitions, due to the 
effective long-ranged interaction of the order parameter. 
We showed that the nodal liquid phase is stable against 
an AF instability within a sizable window of AF exchange 
coupling and that there exists a non-trivial fixed point 
at the AF transition. Assuming that it is possible to 
access the quantum critical regime, critical behavior of 
thermodynamic quantities would be observed provided 
that the distance from the quantum critical point t can 
be controlled. In principle, t depends on the AF ex- 
change interaction which in turn is controlled by dop- 
ing. Estimates from two-magnon Raman scattering ex- 
periments on cuprates [O show that the effective AF 
super-exchange coupling J takes a value of ~ 125 meV 
for AF insulators and decreases as doping. Here we as- 
sume that the only effect of doping is to reduce the effec- 
tive exchange coupling, and also ignore the possibility of 
a spin-density-wave ordering with incommensurate wave- 
vectors Jl5| . We expect a scaling form of thermodynamic 
quantities in the quantum critical regime such as the crit- 
ical specific heat coefficient, j v = c v /T ~ T d ~ 1 F(t vz /T) 
in terms of a scaling function F(x). From the scaling rela- 
tion, we find that zv = .27/(2-77) w l+0.78/n+O(l/n 2 ). 
This analysis may provide a check on the notion of the 
nodal liquid as the zero-temperature pseudogap phase of 
the underdoped cuprates. Also we notice that the d-wave 
like excitation is possibly ubiquitous even in undoped or 
underdoped insulating cuprates, as recently revealed by 
the angle-resolved photoemission spectroscopy on an in- 
sulating parent compound Ca2Cu02Cl2 and theo- 
ries on nodal excitations merit further attention. 

It would be also interesting to study the effect of nodal 
quasiparticles on the superconductor-insulator transition 
in the cuprates. At first glance, however, the coupling of 
the superconducting order parameter to nodal quasipar- 
ticles is suppressed by a d-wave angular factor of cos2</> 



near <f> = 7t(2to + 1)/4 and we expect that it is subleading 
in the RG sense near the quantum phase transition of the 
2 + 1 dimensional XY type. 
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